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551.558.21 
A COMPUTING AID TO STUDIES OF AIRFLOW OVER MOUNTAINS 
By C. E. WALLINGTON* 


Summary. A method is described for using a high-s computer to compute and plot 
flow patterns and orographic cloud over and in lee of given high-ground profiles. Care is 
needed in making boundary assumptions and in interpreting the results, but the method can 
be a powerful aid to observational studies of mountain airflow, including lee waves, rotor 
flow and, possibly, local blocking and high-level turbulence. 


Introduction. In two papérs in 1949 and 1953, Scorer’:* discussed the 
occurrence of lee waves in the atmosphere and showed how the values of 
possible lee wavelengths may be calculated for an airstream that can be 
divided into two or more layers each with a constant value of a parameter 
depending on the wind and temperature structure. The calculations become 
very lengthy if the number of layers exceeds three, so early calculations were 
restricted to simplified airstreams with only a few distinct layers. Corby and 
Wallington’ showed that in some situations lee-wave amplitudes are very 
sensitive to small changes in the airstream, and although simplified models 
may be adequate in many situations there is no simple criterion for assessing 
whether such simplification is justifiable for any particular airstream. There- 
fore, when calculating lee-wave amplitude it is wiser to overcome the computing 
problem than to apply the theoretical technique to a crude approximation 
of a real airstream. 


Use of a high-speed computer for numerical studies of lee waves was 
described by Wallington and Portnall*:® in 1958. This work was exploratory 
and programmes were not designed for general purpose numerical study of 
the lee-wave features of mountain airflow. With a modern computing system 
it is possible to compute and present some mountain airflow calculations in 
such a way that they can be used not only to study lee waves but also as an 
adjunct to studies of rotor flow, local blocking by mountain ridges and high- 
level turbulence. This paper describes the operational features of a computing 
programme that has already been written and indicates why the use of such a 
programme and interpretation of results must be accompanied -by an under- 
standing of the physics and formulation of the problem. 


Formulation of the problem. We shall consider only small perturbations 
in frictionless, steady, isentropic flow in two dimensions in a vertical plane 
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and limit the discussion to waves that are short enough for the Coriolis force 
to be neglected. This is not the only way to formulate the problem, nor is it 
the best way, but it is a method that has already proved to be practicable 
and to produce useful results. If } is the perturbation displacement of a 
streamline from its original level z, then the gas equation, the adiabatic 
equation and the equations of motion and continuity can be manipulated to 
show approximately that 


a*y ‘2 OU 1 #\ dy : Ry ee 

az? t (3 bz P dz dz +GA =) az? = O, in wee 

where U is the undisturbed component of the wind in the x direction at 
height z, x being the distance along the horizontal axis of the vertical 
cross-section, 

® is the undisturbed component of the potential temperature at 
height z, 
and g is gravitational acceleration. 


In order to seek solutions which are sinusoidal in the x direction we 
substitute | = the real part of f(z)e**, where f(z) is a function of z only, 
k represents wave numbers to be determined, and i = »/-1. Equation (1) 
may then be written in the form 


e*f(z) 
Ao 


i> 
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25 ote) + (2 +k) F(z) = 0, ++ (2) 


where § and / are functions of z only. 
If we follow Scorer and consider the flow across a mountain ridge whose 
height, A, is given by 


Hb? 
h= ——. aging 
eng (3) 
where H is the height of the crest and 6 is a width parameter, then the solution 
for the streamline displacement becomes 


k= @ 
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k=0 
where f(o) denotes the value of this function at ground level, z = o. 


The integral will have singularities for any values of k that yield f(0) = o. 
These values of k (= K, say) are the lee-wave numbers, and the flow pattern 
corresponding to them are called the lee waves. On the assumption that 
these lee waves do not occur upstream of the mountain ridge that gets them 
in action, equation (4) may be written 


Y = real part of Hb | e~ k(b + ix} te dk - 
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where the integral is taken around a circuit that excludes singularities, and 
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the second term is repeated for all wave numbers but applies only to the 
downwind side of the mountain ridge. 


The first part of the solution involving the contour integral is called the 
barostromatic flow, while the second part refers to the lee-wave flow. 


It can be shown that if b> 1// at all heights, the integral for the barostromatic 
flow is dominated by the contribution at k =o. Thus the barostromatic 
displacement tz, is given oT eee by 


be = real part of = ——. (b + ix) (4)... , : ee 


This solution includes an oe as the boundary cain can be 
satisfied by either the positive or negative sign before the ix. The effect of 
the ix in the equation is to shift the crests or troughs in the barostromatic 
streamlines over a mountain ridge in the upwind or downwind directions, 
but for the present it will be assumed that the tilt is small enough to be 
neglected. Then equation (6) becomes 


f(z) 
ag a - (7) 

The flow pattern across any high-ground cross-section can be determined 
by analysing the cross-section into a number of ridges of the type specified in 
equation (3), computing the flow across each of these ridges, and adding these 
flows together to form the complete pattern across the high-ground cross- 
section. A computing programme has been written to calculate flow patterns 
in this way, and to compute, as by-products, wave-flow parameters, such as 
the / in equation (2), and the effects of the wave flow on radiosondes ascending 
through the disturbed airstream. 


ve = height of ground x real part of ( 


The computing programme. Data required for the programme 
comprise : 

(i) Airstream data. Winds, temperatures and dew-points at as many levels 
as are available. The levels at which these items are given can be in millibars, 
feet, kilometres, or metres; the programme deals with any of these units. 
Wind directions are in degrees; speeds are in kilometres or miles per hour, 
metres per second, or knots. Temperatures are in degrees Kelvin, Celsius 
or Fahrenheit; dew-points are not essential if cloud computations are not 
required. 

(ii) High-ground data. A list of heights at specified regular intervals on a 
cross-section of the high ground in a specified direction. Alternatively, if an 
experimenter wishes to study the flow over one or more ridges of the type 
specified by equation (3), he need only give the height, width and position 
parameters of the ridges. 

(iii) Boundary levels and number of computing levels. The heights of the upper 
and lower boundaries between which the airflow calculations will be made 
and the number of levels between these boundaries that will be used for 
finite differences in the vertical finite difference form of equation (2). 

(iv) Release points and rates of ascent of simulated radiosondes. If the experimenter 
wants the programme to calculate the temperatures, dew-points and wind 
speeds that would be measured by a radiosonde ascent in the flow, he must 
specify the radiosonde release point and the rate of ascent. 
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With the data in (i) to (iv) the computing programme goes through all 
or a selection of the following actions : 
(a) ‘The airstream data are printed out and plotted in graphical form on a 
plotter. Figure 1 shows a typical plot. 
250° 





\ 
sitatoserr id Bitsriasy 
00 200 kmh 


FIGURE I-——-PROFILES OF TEMPERATURE, DEW-POINT AND WIND COMPONENT 


FROM 270° FOR 0950 EST, 17 APRIL 1966 AT HOBART, TASMANIA 
The shaded layer is the cloud assumed likely if the dew-point is within 1 degC of the temperature. 


(b) Temperatures, dew-points and wind components across the high 
ground for the levels to be used in vertical finite differences are calculated 
by interpolation from the airstream data. 

(c) The parameters S$ and / for equation (2) are computed and graphs 
of / and 2n/l are plotted against height. 

(d) Starting with arbitrary values (to be discussed in the next section 
on page 162) at the top two levels, the simplest finite difference form of 
equation (2) is integrated down to the lower boundary with k = o. 

(e) Starting with similar boundary conditions to those just described, 
lee-wave numbers, A, are found by the method described by Wallington 
and Portnall.* 

({) Graphs of the barostromatic displacement factor and 
an[f/z)/(@f(o)/@k) lx x are plotted against height. 

(g) The high-ground profile is analysed into ridges of the form specified 
by equation (g) by the method described by Wallington. * 

(h) The heights of 20 streamlines between the upper and lower boundaries 
are calculated such that they represent the undisturbed wind components, 
i.e. the wind speed is proportional to the streamline gradient. 
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(t) The following calculations are made for each streamline. For each 
point over the given ground cross-section a search is made from ground 
level upwards to find the height, or heights, at which the streamline should 
be to make its displacement equal to the sum of the barostromatic displace- 
ment and the displacements due to lee waves from all upstream ridges. 
The streamline is drawn on the plotter in a vertical cross-section showing 
the flow over the high ground. 


(j) If the airstream data include dew-points, the programme computes 
and draws the cloud pattern in the airflow on the assumptions that cloud 
forms if the temperature is within 1 degC of the dew-point and that 
descending air.in a layer of cloud does not become unsaturated until it 
descends below the base of the layer. 


(k) If the data include release points and rates of ascent of simulated 
radiosondes, the paths of the sondes are calculated and drawn on the 
cross-section, graphs of the temperatures, dew-points and wind speeds 
that the sonde would measure are plotted, and the vertical speeds of the 
air through which each radiosonde ascends are plotted for each level used 
in the computations, 






































4 ' T T 
Hobart 100 


FIGURE 2-—-COMPUTER-PLOTTED STREAMLINE AND CLOUD PATTERN FOR THE 
AIRSTREAM REPRESENTED BY THE DATA IN FIGURE I FLOWING OVER A HIGH- 


GROUND WEST-TO-EAST CROSS-SECTION THROUGH HOBART 
The clouds were shaded in by hand but their outline was drawn by the computer plotter. 
The computed lee wavelength of 18-5 km (10 n.miles) is in good agreement with a wavelength 
of 11 n. aie observed by satellite photography. 


Figure 2 shows the airflow and cloud pattern computed for a flow whose 
wind and temperature profiles are illustrated in Figure 1. The clouds were 
shaded in by hand after the computation but the cloud outlines were drawn 
by the computer-controlled plotter. Airstream winds and temperatures for 
this computation were those measured by radiosonde at Hobart, Tasmania, 
at og Est (GmMr + 10 h) on 17 April 1966. This airstream produced bars 
of lee-wave clouds clearly visible on a photograph taken by Essa 2 satellite 
on orbit 600 at 0740 Est on 17 April 1966. Andersen® described the occurrence 
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and measured the wavelength from the photograph as 11 nautical miles. 
The computed lee wavelength was 10 nautical miles, which is in good agree- 
ment with the observations. 


The boundary assumptions. The experimenter can choose the upper 
and lower boundary levels for the computed flow pattern. The lower boundary 
is taken to be either ground level or some surface which is taken to be an 
effective lower boundary, e.g. if winds are practically calm in a shallow layer 
of stable air below, say, a few hundred metres, the experimenter may wish 
to try a calculation on the assumption that the lower boundary is virtually 
the top of this almost stagnant air together with whatever high ground 
protrudes up through it. In some computations streamline displacement 
is found to be zero at a level below the top of the high ground in the region. 
In such a situation it must be assumed either that the theory is inapplicable or 
that the high ground will block or divert the low-level flow. If a level of zero 
displacement is below the crests of a mountain range but above the level of 
some of the passes, it may be argued that the airflow will be concentrated as 
particularly strong winds through these passes. But this is conjecture at 
present. The main points being made here are that in some situations the 
choice of a low-level boundary is a subject for experimental and observational 
studies, and that the computing facility is an aid to such studies. 


In setting the upper boundary conditions the experimenter can set the 
height and the value of f(z) at the top two levels. By setting these values of 
f(z) he virtually sets the ratio of the values at the top two levels; the 
magnitudes of the displacement values are scaled by the lower boundary 
factors in the displacement terms. There is no a priort reason why any particular 
set of upper boundary conditions should be correct, and no special justification 
for assuming that displacements will become negligibly small at or towards 
some great height. In this numerical method the procedure is to choose one 
or more sets of upper boundary conditions arbitrarily, and then consider the 
nature of the computed flow patterns. If the computed displacements close 
to the upper boundary turn out to be very small compared with the low- 
level displacements, or if the low-level flow appears to be insensitive to 
variations in the choice of upper boundary conditions, it can be argued that 
the computed flow pattern at low levels is as valid as the formulation of the 
problem. Otherwise it must be accepted that the computed flow pattern 
cannot be justified without justification for the particular upper boundary 
conditions chosen. Wallington and Portnall‘ have discussed the effect of 
the upper boundary assumptions on lee-wave calculations in detail. 


The upper boundary effect on barostromatic displacement is similar to its 
effect on lee-wave computations in that there appear to be many situations 
in which the low flow is not very sensitive to variations in the assumed 
conditions at some upper boundary. But the barostromatic computations 
are more complicated because the imaginary part of f(z) must be considered. 
If the real part of f(z) in equation (6) is almost independent of the assumed 
upper boundary conditions it will have a vertical distribution of shape similar 
to that of the imaginary part. Therefore, the factor [f(z) / f(0)]« — 9 will be 
almost entirely real, in which case the crests or troughs in the barostromatic 
flow will have negligible tilt from the vertical over the mountain ridge and 
the ambiguity of sign is unimportant. 
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However, if the flow pattern is sensitive to variations in upper boundary 
assumptions the real and imaginary parts of f(z) will not necessarily have 
similar profiles and the computed magnitudes and precise locations of troughs 
and crests over or near the mountain ridge will both be of doubtful accuracy. 


Rotors. Because the displacement calculated for a level z refers to displace- 
ment at that level and not from that level it is possible for a streamline to 
have more than one level over any point on the ground profile. Ifa streamline 
has three possible positions it must have either an S shape or be part of a closed 
circulation in the vertical cross-section. It must be stressed that the perturba- 
tion assumptions are not strictly justified where a local reversal of flow occurs. 
However, comparison of computed and observed patterns suggests that a 
local breakdown does not necessarily invalidate the whole computed flow 
pattern. A computed rotor or reversal in the flow pattern can be taken to 
indicate a local breakdown of the flow into chaotic or unsteady flow, rather 
than a simple circulation in a vertical plane. Figure 3 shows a computed 
pattern with reversed flow and rotor features. When this situation occurred 
turbulence below about 2 km was so violent that light-aircraft flying was 
abandoned for the day. Turbulence associated with this type of rotor flow 
is not produced merely by a vertical circulation; it is the result of the smaller- 
scale instability that is generated when overturning of air in the stable layer 
associated with the wave flow forces warm air under cooler air from below. 














——— 


FIGURE 3—COMPUTED FLOW PATTERN OVER A WEST-TO-EAST CROSS-SECTION 
OF THE BRINDABELLA MOUNTAINS, AUSTRALIA 


Wind and temperature data used were those measured by radiosonde at °9 EsT, 12 Jul 
at Wagga Wagga, New South Wales. If streamlines had been drawn at closer Red me oe they 
would have shown additional closed circulations, or rotors, under the two wave crests near 
the centre of the pattern. The airstream contained a shallow temperature inversion of 5 degC 
between goo and 890 mb. 


Figure 4(a) shows a well-developed wave and rotor flow investigated during 
the Sierra Wave Project. This has been described by Holmboe and Klieforth.’ 
Figure 4(6) shows the computed flow pattern for the situation. The airstream 
data used were those from Merced, California, at 13 pst (21 GMT), 16 February 
1952. The computed pattern is in good agreement with the observations; 
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<-— 20 km ———_—__--__-> 
FIGURE 4( a ) —~FLOW PATTERN DEDUCED FROM OBSERVATIONS AND MEASUREMENTS 
ON 16 FEBRUARY 1952 
The observed lee wavelength was 18 km at about the 3-km level. 











FIGURE 4(6)—-COMPUTED FLOW PATTERN FOR 16 FEBRUARY 1952 
Airstream data used were those obtained at Merced at 13 pst (21 out), 16 serra 1952 
rotor circulation, bu’ 


The computed streamlines do not include be-endepin® tnautine.o abet t such 
a circulation beneath the wave crest in more-detailed printed results. The 
orographic cloud pattern is omitted for clarity; it contained wave clouds in the correct locations 


with tops at 5 km, dhe he taterenns teenll avecteer. 
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both observed and computed lee wavelengths were 18 km, and although the 
plotted streamlines do not happen to include a closed isopleth, rotor-flow 
circulation below the wave crests was evident from more-detailed printed 
results. 


Some computed flow patterns have rotor-flow features at high levels. These 
may well be indicators of high-level turbulence; the computing technique can 
be an aid in investigations of such turbulence, which can be favoured rather 
than inhibited by static stability. But in the computations made so far there 
is not enough evidence to assert whether the high-level rotors are likely to 
be real or spurious by-products of upper boundary assumptions. 


Conclusion. Care is needed in making boundary assumptions and 
interpreting the results of a mountain airflow computation, and there are 
occasions when an experimenter must not be tempted to place undue con- 
fidence in a computed flow pattern; but, if this need for caution is realized, 
the perturbation and computing technique described here can be a powerful 
aid to observational studies of mountain airflow including lee waves and 
rotor flow and, possibly, local blocking and high-level turbulence. 
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METEOROLOGICAL NOTES ON TRAVEL OF DEBRIS FROM AN 
ATOMIC EXPLOSION AT MURUROA ON 7 JULY 1968 


By C. L. HAWSON 


Summary. The likely dispersal of debris from the first of a series of atomic explosions at 
Mururoa in July 1968 is considered, including spread by eddy diffusion. Times of first detection 
at Pretoria, Aspendale and Bombay of debris from this explosion are compared with tracks 
and times inferred from charts of average upper winds over the world for July and their 
associated standard vector deviations. Good agreement is found. In particular, rapid inter- 
pany i exchange in the lower troposphere is in accord with recent studies of wind fields 
over 


Indian Ocean. 

Introduction. Atomic explosions took place at Mururoa, 22°S 139°W, at 
2200 GT on 7 July 1968 and on subsequent dates from 15 July to 8 September. 
Debris from these explosions was first observed on detectors exposed at ground 
level at Pretoria, 26°S 29°E on 16-17 July, Aspendale (Melbourne), 38°S 
145°E, on 15-22 July and Bombay, 19°N 73°E on 23 July (Peirson’). 

This note compares these observations with tracks and travel times deduced 
from average winds and their associated standard vector deviations at various 
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levels (Heastie and Stephenson,* Tucker*) and from» papers’ concerning 
winds at low levels over the western Indian Ocean (Findlater**). 


To avoid ambiguity the term ‘average wind’ is used to indicate an average 
over a long period of time of the winds at a level, and the term ‘mean wind’ 
to indicate the mean through a layer. Hence, ‘average mean wind’ denotes 
the average over a long period of time of the mean wind in a given height 
interval for a given track. 


Basic meteorological factors affecting the spread of debris. Debris 
from atomic explosions of sufficient size can be expected to rise quickly to 
near the top of the troposphere (about 16 km at 22°S) at least, and then 
move with the winds encountered at each level. Subsequently debris moving 
almost horizontally at any one level can also change level by various mechan- 
isms, including gravitational settling, turbulent or eddy diffusion, and vertical 
motion of the air in organized weather systems. 


In a cloud formed by an atomic explosion a very wide range of particle 
sizes would be generated and in the circumstances of the explosions under 
consideration it is reasonable to assume that a plentiful supply of a wide 
range of particle sizes was injected into the high troposphere (16 km), sufficient 
to take all times between hours and months to fall through the air under 
gravity to the ground. Turbulent or eddy diffusion, like molecular diffusion, 
acts to transfer matter from volumes of higher average concentrations to 
volumes of lower average concentrations. The diffusion coefficient which 
can be used to describe this transport when meteorological effects are included 
is very much larger than the molecular one, and the horizontal diffusion 
coefficient is very much larger than the vertical one. Values of latitudinal 
horizontal eddy diffusion coefficients, K, of about 10'° cm?/s seem likely in 
the free atmosphere in tropical regions, and values of about 10° cm/*s are 
likely for the vertical eddy diffusion coefficient in the troposphere (Bolin’). 
With these values the depth (kilometres) and width (degrees of. latitude) 
of a layer in the atmosphere within which 50 per cent of radioactivity injected 
at a point would be found after various times is shown in Table I. It must be 
appreciated that these eddy diffusion coefficients and associated depths and 
widths are based on long-term average behaviour-and so are ‘only coarse 
estimates for a specific occasion. 


TABLE I-—-DEPTH AND WIDTH OF A LAYER IN THE ATMOSPHERE WITHIN WHICH 
50 PER CENT OF RADIOACTIVITY INJECTED AT A POINT WILL BE FOUND, IF THE 
DIFFUSION COEFFICIENTS ARE ASSUMED TO BE AS SHOWN : 
Days after DEPTH WIDTH 
injection K = 105 cm?/s K = 10" cm?/s 


— oer. | 


are 
7 
10 
I 
16 pS 


32 28 
XK = diffusion coefficient . 1 degree of latitude 2 111 km -- 


Vertical motion arising from organized weather patterns may be upward 
or downward. . Anticyclones are characterized by slowly descending air 
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(subsidence) throughout a large part. of the troposphere. The subsidence 
proceeds at a maximum rate in the early stages of formation of the anticyclone 
and is greatest in mid-troposphere when values of about 1 km/day may be 
reached, but they are usually very much less. Cloud and rain are associated 
with upward vertical air motion and for widespread rain average velocities 
are usually small, about 5 to 10 cm/s in moderate rain; but in thunderstorms 
violent upward and downward air motion occurs with speeds of 5-10 m/s or 
more. Photographs from satellites show that the relevant areas between 
Mururoa and South Africa and Australia were mostly free of cloud during the 
period 7-22 July 1968. This indicates an absence of significant upward 
organized vertical air motion and probably very slow anticyclonic subsidence. 
Any such slow descent, added to the descent due to gravitational settling 
through the air of the very wide range of particle sizes assumed, would change 
the particular particle sizes which would descend from a specific height in a 
specific time and vary the relative times such particles spent in specific layers 
of the atmosphere (increasing the time in the high troposphere relative to 
that in the middle troposphere), but would not affect the assumption that 
particles would take all times between hours and months to fall from the 
high troposphere to the ground. Indeed, the absence of cloud suggests that 
aggregation processes would be operating in unfavourable conditions so that 
long-time fallout from the very small particles would be prolonged beyond 
the average. 


Discussion. Charts of the 5-year average wind fields during July at 
100 mb (16 km), 150 mb (14 km), 200 mb (12 km), 300 mb (9°5 km), 500 mb 
(5°5 km), and 700 mb (3 km) are available.*;* These show Mururoa to be 
situated near the northern edge of a broad belt of westerly average winds 
which extend right round the southern hemisphere at all these levels. Stream- 
lines from 22°S. 139°W pass across Africa between latitudes 30° and 22°S 
for all these levels from 100 to 500 mb, and across western Australia between 
latitudes 18° and 22°S. At 700 mb, streamlines over southern Africa and 
over Australia are westerly with small perturbations. Some sample average 
wind speeds (all the average wind directions are from the west) and associated 
standard vector deviations of wind at the various levels, read from the charts 
for Mururoa, Pretoria and Aspendale, are given in Table II. 


TABLE II-—-AVERAGE WIND SPEEDS AND ASSOCIATED STANDARD VECTOR DEVIA- 
TIONS OF WIND AT STANDARD LEVELS DURING JULY 
Level Mururoa Pretoria Aspendale 


Average s.v.d, Average s.v.d. Average s.v.d. 

mb - knots 

100 16°5 5 30 42 25 
150 14 & bo 37 5° 35 
200 12 60 60 3 50 40 
goo 9°5 55 50 44 38 40 
500 5°7 33 20 25 25 30 
700 3 <10 <10 20 10 22 


Values are extracted from Geophysical Memoirs Nos 103% and 105.3 Average wind 
directions are all westerly. , 
At 200.and 150 mb average winds to the south of Mururoa increase by up to 20 kt 
within 2° or 3° of latitude. 

The mean tropopause pressure in the area of Mururoa is close to 100 mb. 

s.v.d, = s' vector deviation. 
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It is assumed (i) that the settling speeds (including the effects of diffusion) 
of radioactive particles arriving at the ground at Pretoria were uniform with 
height and (ii) that the particles commenced their fall from various heights 
over Mururoa and moved with the average winds for the heights and locations 
encountered en route. It follows that the mean average westerly wind com- 
ponents are 35 kt for the layer 16 to o km, 32 kt for the layer 14 to o km and 
less for layers starting at lower levels and finishing at the ground. Values 
for layers from levels greater than about 18 km to the ground would be less 
than 35 kt. Values for layers to 3 km above the ground would also be a 
maximum from a level of about 16 km, for which 43 kt is indicated. 


If it is assumed that the filter at Pretoria was placed in position during the 
morning of the 16th and removed 24 hours later, the elapsed time between 
the first explosion at Mururoa and the first detection at Pretoria is between 
84 and 9} days. The maximum mean westerly wind component Mururoa— 
Pretoria indicated by the radioactive debris transit is thus 45-41 kt, corres- 
ponding to travel times of 8} to 9} days respectively. Sample values of 
standard vector deviation of the winds at the standard levels during July 
are given in Table II. The standard vector deviation of mean wind cannot 
be inferred directly from the appropriate values along tracks since it depends 
also on the correlation of wind with height, distance and time. These are 
not specifically known but, considering the manner in which correlations 
change with these three factors, it is plausible to assume in the circumstances 
here that although winds in narrow layers will be highly correlated the 
correlation will fall off very rapidly with increasing layer thickness (Durst,* 
Kochanski*®). Assumptions of correlations of 1 over layers 2 km thick and of 
o between such layers suggest a standard vector deviation of mean winds 
16-o km for the track Mururoa—Pretoria of about 10 kt. With the further 
assumption of a circular normal vector distribution about the average mean 
wind of 35 kt, a 13 per cent chance of encountering a mean westerly wind 
component of 43 kt or more is indicated for the 16-0 km layer Mururoa-— 
Pretoria. The value obtained from the radioactive debris transit Mururoa— 
Pretoria thus readily accords with the average charts and the assumptions 
made that the debris involved started from the top of the troposphere and 
fell steadily (in the mean) to the ground. An alternative explanation that the 
debris fell more slowly in the high troposphere than elsewhere is also possible. 


Considerations similar to those discussed above, applied to the track 
Mururoa—Aspendale, yield an average mean wind 16-0 km of 37 kt, a standard 
vector deviation of about 1o kt, and probabilities of 18 and 4 per cent that 
the mean 16-0 km wind exceeds 44 and 51 kt respectively. Radioactivity 
first appeared at Aspendale upon a filter exposed from 15 to 22 July; this 
indicated minimum travel times of 7 to 14 days corresponding to mean westerly 
wind speeds of 88 to 44 kt respectively (51 kt for 12 days’ travel). These 
figures suggest that the radioactive debris first arrived at Aspendale towards 
the end of the period covered by the filter, began its voyage near the top 
of the troposphere and fell for the most part steadily to the ground. The 
departure of the latitude of Aspendale from that of the mean latitude given 
by the average streamlines is some 18° and appears rather large. However, 
the concentrations of barium-140 observed at Aspendale are appreciably 
lower than those observed at Pretoria and the ratio of the two concentrations, 
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appropriate to the times of first detection at the two stations, accords quite 
well with those to be expected from considerations of eddy diffusion about 
the average mean tracks 16 km to the ground which were assumed. 


Turning attention to the first arrival of debris at Bombay it is convenient 
to begin by considering back tracks from Bombay using Findlater’s average 
charts below 3 km. These indicate a curving track over Somaliland and 
thence south and south-eastward to an area east of Madagascar and average 
mean wind speeds between about 20 and 30 kt, corresponding to travel times 
to the region of 15°S 60°E of 9 and 6 days respectively. By comparison 
between the average mean winds 16-0 km and 16-3 km obtained for the track 
from Mururoa to Pretoria, the mean 16-3 km wind from Mururoa to 60°E is 
likely to be 5 to 10 kt greater than the mean of 43 kt indicated by the debris 
travel time from Mururoa to Pretoria. This suggests that the 3-km level around 
15°S 60°E probably first became polluted 9 days after the explosion. This 
pollution gradually falling to lower and lower levels below 3 km would then 
be likely to follow the track explored above to reach the ground at Bombay 
some 6 to g days later, giving a total minimum travel time to Bombay of 
15-18 days. This total is in good agreement with the observed first appearance 
at Bombay on a filter exposed on 23 July, about 15} days after the explosion. 
Whilst there seems little doubt that this route to Bombay is probably close 
to that actually followed by the debris, it should be noted that the average 
descent rates used in the two parts of the argument are different, i.e. descent 
from 16 to 3 km in g days (1-4 km/day) and from 3 to o km in 6 days (0-5 
km/day). This can be qualitatively justified by the assumption of subsiding 
air in anticyclonic conditions in mid-troposphere, as discussed earlier (p. 166). 
Debris from the explosion following substantially the suggested route to Bombay 
and passing on eastward in the lower troposphere would be expected to be 
largely washed out in the Indian monsoon. Other trans-equatorial flows 
exist at relatively low levels that would feed debris into the northern hemisphere 
from the broad belt established around the world in the tropical regions 
of the southern hemisphere. However, these would not usually be expected 
to carry air far from the equator before pollutants were washed out by rain 
associated with these currents. The arrival of debris at stations in the northern 
hemisphere remote from Bombay indicates that debris was carried across 
the equator by flows other than the Findlater current during this period. 
Unfortunately investigation is complicated by the succession of explosions 
and the broad periods for the times of arrival of the debris at the stations, 
so that detailed study becomes unprofitable. 


Conclusions. The times of first detection of radioactive debris at the 
ground at Pretoria, Aspendale and Bombay following an atomic explosion 
at Mururoa on 7 July 1968 accord well with times and tracks inferred from 
charts of average upper winds over the world for July and their associated 
standard vector deviations. In particular, rapid interhemispheric interchange 
in the lower troposphere over the Indian Ocean is in agreement with recent 
studies by Findlater of the wind fields in that area. 


ents. The author wishes to thank Dr Peirson of the Atomic 
Energy Research Establishment, Harwell, who brought the observations to 
his attention, and Mr R. F. Jones, Meteorological Office, Bracknell, for helpful 
advice and encouragement. 
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AGROMETEOROLOGY IN THE METEOROLOGICAL OFFICE 


By G. W. HURST 


Summary. A survey is given of the work and ization of the Agricultural Branch of the 
Meteorological Office over the last five & ). Reference is made to the scope of 


research and field ex ts un also to several co-operative investigations 
(mostly with the ey a. ture, Fisheries and Food) into major agricultural and 
horticultural probl tion is given of future trends. 


Introduction. A short account was given some time ago (Hurst’) of 
services which the Agricultural Branch gave to the agricultural community, 
based on a colloquium on the subject in November 1964. In five years much 
can happen, and it is perhaps appropriate to review the present situation, 
particularly in relation to research projects which have been, or are being, 
conducted. 

The Branch still continues to work closely with the Ministry of Agriculture, 
Fisheries and Food (MAFF) in England and Wales and, in particular, with 
its National Agricultural Advisory Service (NAAS), which has eight Regional 
Headquarters including those at Bristol, Cambridge and Leeds at each of 
which two or three meteorological staff are based. Of the others, Reading 
is supported from Bracknell, and Aberystwyth, Derby, Newcastle and 
Wolverhampton at present rely for advice on the meteorological staff at 
Bristol, Cambridge or Leeds as geographically appropriate; but it is hoped 
that eventually agrometeorologists will be located at all eight Headquarters. 
The task of such a meteorologist is to advise the Regional Director, to assist 
in solving problems raised by his staff, and to consider the application of 
meteorology to agricultural and allied problems. This type of organization 
does not obtain in Scotland and Northern Ireland. Inquiries in Scotland 
are handled by specialists at the. agricultural colleges, who, as.:necessary, 
consult the agrometeorologists in Edinburgh. Inquiries in Northern Ireland 
are dealt with by the Meteorological Office, Belfast, with some support from 
Edinburgh. 
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Close contact: is also maintained, from. both Bracknell and outstations, 
with the various agricultural research organizations — some of which come 
under MAFF, and some under the Agricultural Research Council (ARC), 
Forestry Commission, etc. — and, of course, with several university depart- 
ments of agriculture and horticulture. In this context, a Senior Scientific 
Officer is:seconded to the Grassland Research Institute at Hurley, Berkshire, 
for a minimum period of two years to assist in problems of meteorology and 
grass production. It is hoped that this will be the first of a series of second- 
ments of this type, probably for periods longer than two years in order to 
allow sufficient time for the development and completion of a worthwhile 
programme. The strength of the Branch numbers under 20 — most of whom 
own gum boots ! 


Services. These were mentioned in the earlier review as falling into two 
classes: routine and non-routine. Little change in routine has taken place, 
e.g. weekly weather summaries are still prepared and sent to NAAS and 
other recipients, and warnings are still issued for potato blight and apple 
scab, though such warnings are increasingly being channelled through NAAS 
(either Regional Centres or the Pathology Laboratory at Harpenden). 
Similarly, weather information related to potential severity of liver fluke 
and other diseases is passed to the Central Veterinary Laboratory at Wey- 
bridge. Transpiration figures are issued for irrigation need. 

Non-routine services include the answering of queries received either direct 
from farmers and others or, more frequently, through NAAS and similar 
channels. Lecture invitations unfortunately have usually to be refused 
because staff does not exist for the purpose, but certain commitments are 
accepted, such as the annual lectures to Wisley, or Kew, horticultural students, 
and the bi-annual symposium for the benefit of the agrometeorological 
observers course at the Meteorological Office Training School. In 1969, 
the first Voluntary Service Overseas Briefing Course was held at Reading 
University and the Branch participated actively in the programme. Again 
unfortunately is is seldom possible to spare personnel to support agricultural 
shows but occasionally an exhibition is mounted (see Plate I); in addition 
to instruments familiar to most meteorologists, this display included a wetness 
recorder, where dew deposited on the ‘toffee apple’ polystyrene knob is 
counterbalanced and a record is made on a barograph type of chart. These 
instruments are’ particularly used by apple and potato growers in spring 
and early summer. 


“Research over the last five years. The area covered has been very 
wide and an article of this length cannot cover all aspects. ‘A representative 
selection is therefore discussed below, with brief indication of what has been 
done and of future. developments, if appropriate. This is perhaps the place 
to mention that during the last 5 years 10 workers in the Branch published 
over 80 papers amongst them. 


- Foot-and-mouth disease. Following: the unprecedented peer of «Hin ta 
ment of the disease after the outbreak in October 1967, the Branch was 
consulted early in 1968 on what meteorological factors might be involved in 
the spread of the disease. Work was soon completed- by one-member of the 
staff on the possibility of spread from the continent; it was shown. that ‘the 
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disease could not have so originated in 1967, though almost certainly it had 
in some earlier years. The main effort concerned the spread of the disease 
within the country, and two members of the staff worked very closely with a 
pathologist at the Central Veterinary Laboratory, Weybridge, and found 
striking relationships between the spread of the foot-and-mouth virus and 
the wind and rainfall during the incubation period. Thus Figure 1 shows 
the spread of an earlier outbreak in Hampshire in January 1967; there were 
seven primary sources and although from 3rd onwards the wind varied widely, 
the sector in which wind was associated with rain was confined to well under 
180°, and subsequent outbreaks were all within this sector. Papers elaborating 
the results of these researches have been published in several journals. 
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FIGURE I—SPREAD OF FOOT-AND-MOUTH DISEASE IN HAMPSHIRE IN JANUARY 1967 


Sources © Infections from 8-22 January 
Figures of sues are Néthiinal Grid refetences, not dagress: 


Evidence was given by three Branch members to the Northumberland 
Committee of Inquiry, an outcome of which has been the creation of a scheme 
for putting a team (including meteorologists) on the spot without delay in 
future outbreaks. 


Irrigation and water balance. Two particularly important events have 
happened; first, the new average transpiration figures were completed for 
the revised MAFF Technical Bulletin No. 16,* which was issued in 1967, and 
second, in 1967 MAFF published an atlas of long-term irrigation needs for 
England and Wales* compiled by a member of the Branch. This atlas 
indicates the irrigation required over varying periods of the summer half 
of the year to ensure that various planned soil moisture deficits (SMDs) are 
not exceeded, considering the needs of the driest year in 20 and also those 
of the 5th driest year in 20. Figure 2, based on a plate from this atlas, shows 
the irrigation required in the 5th driest year in 20 to restore the ground to 
field capacity every time a 2-inch SMD is reached in the 4-month period 
May to August. 
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PLATE ll(a)—IMPACTOR USED FOR SPORE COLLECTION BY AIRCRAFT 


The right half of the instrument projects horizontally out of the aircraft, perpendicular to the 
side, in free airflow. See page 175. 


Photographs by courtesy of Royal Aircraft Establishment, Farnborough 


PLATE 11(b)—SAMPLING HEAD OF IMPACTOR 


The orifice on the side of the head is protected by a cover (lower left) before operation. Spores 
are impacted on sticky transparent tape wrapped round the drum (centre). See page 175. 





To face page 173 


(a) Whole site of 17 beds. 


Photographs by J. Cochrane 


(6) Single bed showing framework for polythene cover. 


PLATE IV—-WATERCRESS BEDS AT FOBDOWN, HAMPSHIRE 


See page 176. 
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FIGURE 2—IRRIGATION NEED IN 5th DRIEST YEAR IN 20 FOR THE PERIOD MAY 
To AuGusT (after Hogg*) 


Shows irrigation required to restore the ground to field capacity whenever a 2-inch soil 
moisture deficit is reached. 


Work continued on rainfall effects on agriculture and included studies, 
principally through SMD considerations, of excess rainfall in autumn and 
winter, with consequent leaching and poaching problems, and of return to 
field capacity in autumn with subsequent flooding, etc. 

Work being done on machinery ‘work days’ is related to water-balance 
considerations; very broadly the approach is to assess the wet days in spring 
and the water balance in autumn. 
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Cropping and meteorological factors. Some work has already been done 
relating cropping to current and preceding weather and this continues to 
be an important and active field of research. An approach has been made 
to consider effective transpiration, defined as the calculated actual trans- 
piration occurring when the SMD is less than 2 inches in the root zone, and 
in particular the problem of effective transpiration and the growth of grass 
has been studied. Some work has been done on the harvesting and quality 
of crops such as hops, wheat and barley; with hop yields for example, the 
main meteorological variables considered are effective transpiration in spring 
and early summer and sunshine in late summer and the beginning of autumn. 
Work is being done relating hop quality to meteorological factors, especially 
temperature, for a period before harvesting ; ‘Plate II shows maximum and 
minimum thermometers in cocoa tin mounts in hop-field trials in Kent. 
A study is just starting to help to assess the suitability of the climate in southern 
England for maize. production. 


An important application of the research on grass growth has led to good 
forecasts being made in March and June of the probable milk yields up to 
the following spring. The most important meteorological factors are the 
soil temperature excess above 43°F (in March) and the effect of June rainfall 
on subsequent hay quality — dry weather at this stage goes with good quality 
hay. Figure 3 shows well the accuracy of the forecast of annual milk yields 
made at the end of June for the year ending in the following March; a 
regression equation includes factors for the milk yield in the preceding months 
April to June, and a term for the June rainfall. 
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FIGURE 3—FORECASTS OF ANNUAL . MILK YIELDS, MADE AT THE END OF JUNE 
FOR THE YEAR _ENDING THE FOLLOWING MARCH 
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Also related to quantity and quality of-creps is the study of honey production 
and meteorology. No short-term forecast. technique has..been found here. as 
the main weather relation is between - ‘high honey yields and-hot weather in 
July (especially) and- August, but examination of temperatures over the. last 
century suggests that future prospects: for hetser: —— —— are. ae wees 
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An approach has also been made to evaluate crop potential by forecasting 
the level of a particular disease or pest. Thus high February temperatures 
lead to a high level of yellows virus disease in sugar-beet.. Rather similarly, 
consideration was given to the activity of the shallot aphid; which can be 


very damaging to strawberry crops; again, high Tateveny pemapaqneunes are 
associated with high pest incidence. 


Animal diseases. Working closely in conjiinction with a véry energetic 
and interested pathologist at the Central Veterinary Laboratory, a member 
of the Branch has found that the incidence of several diseases’ of sheep and 
cattle is closely related to meteorological factors, e.g. liver fluke in sheep to 
summer rainfall. Examples include pregnancy toxaemia of éwes, the incidence 
of which was found to be higher in cold wet winters; a numerical relationship 
with weather was established for this as for the other diseases. Nematodirus of 
sheep was found to be closely related to the 1-foot éarth temperature, and 
the hatching out of the larvae was rélated to March temperatures, ‘high 
temperatures being associated with an early hatch. A useful p 
forecast was possible in thé second half of March. Gastro-enterjtis of cattle 
is associated with late summer or early autumn rainfall and this, with know- 
ledge of the previous year’s disease figures, enables a fairly accurate assessment 
of the following year’s disease level to be made. 

Other diseases were tnvertignteds such -as pumyieacle of ewes and stillbirth 
in autumn calvers. 


Airborne movement. Much ‘has been done in the examination scape movements 


of insects, spores, viruses, pollens, pollutants, etc.; reference has already 
been made to the movement of foot-and-mouth virus. Many analyses have 
been made using synoptic charts (usually in collaboration with entomologists 
from Rothamsted Experimental Station, Harpenden, or elsewhere), for 
example analyses of particular flights across the Atlantic, and definite patterns 
of immigration have been established for some species of insects. In this 
connection, a monthly analysis of trajectories to Great Britain from the 
continent has been prepared from 1946 onwards using a weather classification 
approach. Analysis has been made and published of spore trapping. over the 
North Sea and elsewhere by the Meteorological Research Flight in conjunction 
with Rothamsted, This has shown very interesting .patterns of spore and 
pollen concentrations in relation to the diurnal temperature over the source 
area (U.K. usually), often with marked maximum counts well away from the 
English coast; vertical changes were considered in relation to stability and 
terminal velocities. The impactor trap, as used on a Hastings aircraft of 
the Meteorological Research Flight, and the sampling head ate shown on 
Plate III. Diffusive processes have also been studied both with the spread of 
potato blight, fireblight, etc. and (using radioactive tracers) with. indoor, 
controlled movement in animal: houses;. movement. in vegetable stores. has 
also been examined. It has been possible to point out SRP wrong 
design of ventilation systems on at least one occasion. : 


Field experiments. The Branch has participated in several important -field 
trials and is still doing so. Three of interest are those in Thetford Chase, 
the orchard trials and the watercress trials. The Thetford experiments were 
made in conjunction with the Forestry. Commission and were. concerned 
with frost damage to young and sensitive Corsican pine saplings. A network 
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of minimum thermometers exposed in ‘cocoa tins’ was maintained for three 
years with modifications in the light of experience; results are proving very 
helpful to the Forestry Commission. The orchard trials were started at the 
request of NAAS to investigate the weather effects on orchard performance in 
different parts of southern England and the west Midlands; the experiments 
were conducted over three years, and preliminary survey of results seems to 
confirm the meteorologists’ expectations that macroscale differences are 
important, but that mesoscale differences seem to matter much less than 
factors such as pollination and husbandry. The watercress trials, started late 
in 1969, also in conjunction with NAAS, are expected to continue for several 
winters. Seventeen similar parallel concrete beds have been laid down at 
a cress farm and temperatures are being measured for different water flow 
rates and for other varying régimes; very interesting results are hoped for, 
linking varied flow rates, cover and crops with fall of temperature as the (well) 
water proceeds down the beds. The general layout of the site at Fobdown, 
Hampshire, is shown in Plates IV (a) and (6). Plate IV(a) shows the whole 
site of 17 beds and Plate IV(6) shows a single bed with cress growing and a 
framework for the polythene cover; probes for mercury-in-steel thermographs 
are located about 2 feet from the ends of several of the beds. 


Land use. Work under this heading mostly does not fall under the title 
of research, apart from certain aspects such as the effect of slope on classifica- 
tion of land grading and the shelter effect of land forms and topographical 
factors generally. A classification is being done by the Agricultural Land 
Services with which the Branch co-operates; it is mainly a map exercise 


with occasional support by field examination. Eventually the whole of 
England and Wales will be covered, using as the basic map the Ordnance 
Survey 1 inch to 1 mile. 

Work in much more detail is also being done in certain areas for a 
horticultural survey of land potential. 


Miscellaneous. Valuable work continues in various directions, but in a 
short space, passing mention only can be made of a selection as follows: the 
relationship of mulching with soil temperatures; the effect of shelter, both 
natural and artificial; meteorological problems in upland climates, including 
continued interest in flag tattering as a measure of wind flow. The Branch 
continues to be represented on the Shelter Research Committee. Important 
pioneering work is also being done on radiation levels within glasshouses, 
and a computer programme is being developed for the assessment of glass- 
house design and location. 


International services. Finally a little must be said of the contribution 
of the Meteorological Office outside the realms of domestic agriculture and 
horticulture. Considering first the World Meteorological Organization 
(WMO) the Branch contribution is distinguished; in addition to providing 
the President of the Commission for Agricultural Meteorology since 1962 
it has provided chairmen for two Working Groups and five rapporteurs on 
various subjects. Six WMO Technical Notes written within the Branch are 
either already published or well under way. 

The Branch has been well represented at conferences. In the five years 
covered, one or more members have attended at least 12 conferences or 
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meetings overseas and many more in this country, and at most of the con- 
ferences, leading papers were presented and sessions and working groups 
were chaired. Pressure of work and other reasons have necessitated refusal 
of some quite pressing invitations from abroad. 


In 1970, plans are in hand for representation at a symposium on ‘Plant 
Responses to Climatic Factors’ in Uppsala (Sweden); it is expected that a 
three-month visit to Turkey will be made by a member of the Branch in 
autumn to complete a six-month period of technical assistance in agro- 
meteorology which started in 1967. It is hoped that in November a Branch 
member will conduct a WMO seminar in Barbados on agrometeorology. 
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NOTE ON THE FORECASTS OF WAVE HEIGHT MADE BY THE 
METEOROLOGICAL OFFICE 


By R. F. ZOBEL, O.B.E. and R. DIXON 


Summary. The method of forecasting wave heights in the North Atlantic Ocean is discussed 
and some results presented. The errors obtained are shown to be very similar to those obtained 
by Pore and Richardson. 


Introduction. The Central Forecasting Office (CFO) of the Meteorological 
Office is designated in World Weather Watch terms, both as the National 
Meteorological Centre and as a Regional Meteorological Centre. Sea wave 
analyses and forecasts are required under both headings, firstly in connection 
with the ships’ routeing service provided by the Meteorological Office and 
secondly as part of the commitments undertaken by a Regional Meteorological 
Centre. Analyses and forecasts of wave height have therefore been produced 
twice a day for periods up to 48 hours ahead and broadcast by radio-facsimile 
since 1 October 1968. An account is given below of the basis of these fore- 
casts and some assessment of the accuracies obtained. 


Derivation of working equation. The requirement is to provide analyses 
and forecasts of the state of sea. This is usually represented by two terms 
consisting of the wind-formed wave and the swell wave. The former is 
dependent on the fetch and time for which the wind has been blowing, whilst 
the latter is related, in a rather complex manner, to events at a distance in 
the recent past. 

Predictions of the height of the wind-formed wave must therefore be based 
on predictions of wind. Swell forecasts should equally clearly be based on a 
series of wind analyses. 


The wind predictions are made in the Meteorological Office on its KDF 9 
computer using an equation derived from data given by Findlater’ et alii. 
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These data may be used to construct curves showing the relationship between 
the goo-mb wind and the surface wind for five ranges of instability in the 
lower layers of the atmosphere. The curves may be converted to curves 
showing the relationship between goo-mb wind and the wave height in feet 
by means of an equation due to Scott. The two curves representative of 
markedly unstable and stable conditions are shown in Figure 1 together 
with a curve representative of the mean of these conditions. The mean lapse 
rate is taken to be about 1-8 degC/1000 feet (1-8 degC/300 m), and it is 
further assumed that the goo-mb wind is adequately represented by the 
geostrophic value obtained frem the computer analyses and forecasts for 
the 1000-mb surface. 
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FIGURE !I-——RELATIONSHIP BETWEEN goo-mb WIND AND WAVE HEIGHT FOR 
MARKEDLY UNSTABLE, MARKEDLY STABLE AND MEAN CONDITIONS 


The actual equations used are as follows : 
- = 10055 - o-0101V,,, + 0°000075V",,,, os» Gt) 


H = 0075 V,*/7 +5, ... (2) 


where H is the height of the wind-formed significant wave in feet, V, is. the 
surface wind and V,,, is the wind at goo mb in knots. The significant wave 
height is defined as the mean height of the biggest one-third of the fully 
developed waves. Full development is yen considered to occur in fetches 
in excess of 400 nautical miles. 


Combination of equations (1) and (2) leads to the ity equation 


H = 0°075 Vogo"/* (10055 — 0°0101 Vg44 + 0-000075 V*4,5)°/? + 5 - «+ (3) 


The constant 5 in equation (2) merely represents the fact that observations 
have shown the average height of the swell wave between 50°N and 60°N 
in the eastern Atlantic to be about 5 feet. This is the only allowance made 
for swell in the system adopted. 


It is of course quite possible to use much more sophisticated techniques 
for both the wind-formed and the swell wave derivations from the existing 
and. pre-existing wind fields, but the philosophy has been that the accuracy 
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of description, especially when predicted, of the wind field is unlikely to be 
sufficiently adequate to justify such additional refinements at least for some 
years. It will appear from the results and discussion below that this philosophy 
is probably sound. 


Findlater et alii‘ found that the angle between the goo-mb wind and the 
surface wind direction averaged about 10°, but on the basis of the philosophy 
mentioned above it was decided to ignore such a value, so that directions 
are taken as those of the computed geostrophic wind. 


Preparation of wave analyses and forecast charts for issue. The 
analyses and forecasts received from the computer refer to wave heights 
which conform to equation (3) and they occur at a regular series of grid 
points irrespective of fetch and, indeed, irrespective of whether the grid point 
refers to land or sea. The forecaster has therefore to make subjective adjust- 
ments to allow for these features and he also makes such other adjustments 
as he feels desirable. For example, the analyses have so far been made without 
any regard to observed wave values, but have been deduced from geostrophic 
wind values with an extremely coarse adjustment for swell. The forecaster 
now pays careful attention to the observations, particularly from ocean 
weather ships, in making his adjustments to the analyses and he also remembers 
that adjustments may be required to allow for cyclostrophic wind effects and 
for marked departures in stability from ‘the average ‘condition assumed in 
developing the working equation. -Figures -2~6 show a computer analysis, 
the forecasts for 24 and 48 hours ahead and the analyses for these times. 























FIGURE 2—-COMPUTER WAVE-HEIGHT ANALYSIS, 00 GMT, 3. JANUARY. 1970 


Wave heights in- metres. 
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FIGURE 3—24-HOUR FORECAST OF WAVE HEIGHT FOR 00 GMT, 4 JANUARY 1970 
Wave heights in metres. 








eee 





FIGURE 4—-COMPUTER WAVE-HEIGHT ANALYSIS, 00 GMT, 4 JANUARY 1970 
Wave heights in metres. 
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FIGURE 5—48-HOUR FORECAST OF WAVE HEIGHT FOR 00 GMT, 5 JANUARY 1970 
Wave heights in metres. 
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FIGURE 6—COMPUTER WAVE-HEIGHT ANALYSIS, 00 GMT, 5, JANUARY 1970 
Wave heights in metres. 
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Results obtained. Table I shows the monthly mean root-mean-square 
errors of the twice-daily computed analyses and 24- and 48-hour forecasts, 
together with the corresponding errors in the charts actually drawn by the 
forecaster. For comparison, the errors in persistence forecasts at 24 and 
48 hours are also given along with the actually observed monthly mean 
combined sea and swell wave. Unfortunately the ‘forecaster’s’ errors are 
only available from August 1969, but it is nevertheless clear that the fore- 
caster’s analyses are more in conformity with the observations than the computer 
analyses are. On the other hand, the errors in the hand-analysed charts might 
be expected to be zero but are not, no doubt because observations of wave 
height are notoriously difficult to make and discrepancies of the kind noted 
are inevitable. It is also clear that the forecaster improves the forecasts 
at both 24 and 48 hours— more at the former than the latter. Both the 
computer and the forecaster beat persistence quite handsomely. 


TABLE I-—-MONTHLY VALUES OF ERRORS IN WAVE ANALYSES AND FORECASTS 
FOR OCEAN WEATHER STATIONS I AND J, DECEMBER 1968 — NOVEMBER 1969 


Ocean Month ROOT-MEAN-SQUARE ERRORS Mean observed 

weather Computer Forecaster Persistence combined sea 

station Analysis Forecast Analysis Forecast and swell wave 
4- -h 24-h 48-h 24-h 48-h 


tw 


feet 
December 1968 
January 1969 
February 
March 
April 
May 
June 
July 
August 
September 
October 
November 


December 1968 
January 1969 
February 
March 
April 

ay 
June 
July 
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MAWUE=— GURMUOwW SCrud6SSs VU rS)—w 
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ABAOUMBA GIABAS BAIGUVHSO OUUWaAY 


VEOGHS SHEDIY SASSAFS GHYUED 
wWOSOAD wWaOuUed 


SSOOYD OHSGAY HG lOO 


Table II shows the same quantities, but as 12-monthly averages, for the 
ocean weather stations A- E and I—K. It will be seen that the results are, 
in general, similar at all stations. The published analyses are much better 
fits to the observations than those from the computer, the forecaster improves 


TABLE II-—12-MONTHLY AVERAGE ERRORS IN WAVE ANALYSES AND FORECASTS 
FOR EIGHT OCEAN WEATHER STATIONS ON THE NORTH ATLANTIC, DECEMBER 


1968 — NOVEMBER 1969 


Ocean weather ROOT-MEAN-SQUARE ERRORS Mean observed 
station Compute Forecaster* i combined sea 
Analysis Forecast Analysis and swell wave 
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on the computer and both are better than persistence. Station A is however 
an exception, in that the computer’s forecasts are better than the forecaster’s. 
The reason for this is not immediately apparent, though it is noted that 
average waves at this station are the highest. 


It is of interest to compare these results with other published values. The 
only figures of a similar kind known to the writers are due to Pore and 
Richardson.* Comparisons cannot be strict as Pore and Richardson only 
give values for May and June 1966 for ocean weather stations A- E, I-J 
and M in the North Atlantic and N, P and V in the North Pacific. Figure 7 
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FIGURE 7——ROOT-MEAN-SQUARE ERRORS OF WAVE-HEIGHT ANALYSIS AND FORE- 
CASTS IN RELATION TO FORECAST PERIOD 


Met. Office computer May-June 1969 Met. Office forecaster 12-month average 
NMC (National Meteorological Center, U.S.A.) May-June 1966 

FNWF (Fleet Numerical Weather Facility, U.S.A.) May-June 1966 

NMC, FNWF and persistence figures are due to Pore and Richardson. 


shows their values and on these have been superimposed the Meteorological 
Office computed results for May and June 1969 together with the forecaster’s 
12-month average values (Table II). It would appear that the values presented 
here and by the U.S. authors are very similar and this may indicate that the 
simple approach to the conversion between forecast wind and wave height 
is justifiable, since the conversion used in the Meteorological Office is the 
simplest used by the three institutions concerned. 
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551-593-553 
NOCTILUCENT CLOUDS OVER WESTERN EUROPE DURING 1969 
By J. PATON 


Accounts of displays of notilucent clouds (NLC) observed over western 
Europe during 1969 are contained in Table I. These clouds are normally 
seen from places between 50°N and 60°N within the period from the end of 
May to early August (though they may very occasionally be visible outside 
this period), so included in the table are notes on observing conditions on 
each night from 31 May to 5 August. 


On nights when extensive tropospheric clouds over western Europe prevented 
the detection of any NLC that may have been present, ‘cloudy’ is entered 
in the notes. When the sky at many stations was sufficiently clear to permit 
the decision that NLC were absent, ‘no NLC’ is entered. 


When NLC were present, the period of duration of reported observations 
of the display is given in the second column, while the third column contains 
particulars of the display including forms and brightness. On occasions when 
it is likely that the cloud field had been observed to be illuminated to its 
southern border at some time during the night, the approximate latitude of 
the border is given.' In the last four columns are given the maximum elevation 
above the northern horizon and the limiting azimuths of the visible cloud 
field recorded at selected stations at stated times during the display. 


TABLE I-——DISPLAYS OF NOCTILUCENT CLOUDS OVER WESTERN EUROPE DURING 1969 


Date— Station Max. Limiting 

night of Times Notes position Time elev. azimuths 
UT UT degrees 

31 May- 


1 June No NLC 
1- 3 June No NLC 
No NLC 
No NLC 
No NLC 
2315-0030 Faint greenish-white bands visible 56:5°N 3°W 
from central Scotland. 


Cloudy 
2102-2345 an ese cary ar oo rk 7 ce 
‘open in part of night. 51° “ 
Later one band predomi: nant afte 
2220 h seen from the 
No NLC 
2325-0140 Faint greenish-white veil 56°5°N 3°W 


2230-0045 Faint bands 55-S°N 4-:S°W 


Bacausa 


No NLC 

2315-0105 Faint cluster of bluish bands. 57-S°N 3-S°W 
Coie S7°N 2°W 
No NLC 


No NLC 
Cloudy 


2345-0145 Veil and ie, Southern boundary 56°N 3°W 
about 58° 55-5°N 1-S°W 


yw 
om 


Cloudy 
No NLC 
Cloudy 
Cloudy 


Cloudy . 
0001-0220 Small patch of bright bands and 5$7:5°N 3-5°W 
billows. 


Southern boundary about 57°N. 56°5°N TOW 
55-S°N 45°W 


54-S5°N OW 
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Date — 


night of Times 


UI 
25-26 June 
2330-2345 


26-27 


27-28 2308-0200 


0140-0150 


2310-0150 


2220-2230 
0140-0150 


2200-2245 


0058-0120 
2130-0015 


2250-0234 


0140-0250 


Notes 


No NLC seen from British Isles. 
Faint NLC seen for short period 
from Denmark (Bornholm). 

No NLC 
Moderately bright bands and billows 
with, around 0005 h, small whirls. 


Southern boundary about 57°N. 


Cloudy 

NLC seen through low cloud from 
Denmark. 

Moderately bright display of bands, 
billows and whirls seen from British 
Isles and over the Altantic. 


Small patch of NLC seen through 
low cloud. 

Cloudy 

No NLC 

No NLC 

Faint NLC seen low 
horizon from Prestwick. 
No NLC 

No NLC 

Faint patches of billows seen through 
low cloud. 


on NNE 


NLC seen from Copenhagen. 

Small bluish-white patch of NLC 
visible only for short period from 
British Isles. 

Cloudy 

Small patch of NLC seen through 
low cloud. 

Cloudy 

Cloudy over British Isles. 

Bands seen close to horizon from 
Denmark (Copenhagen). 

No NLC scen from British Isles. 
Greenish parallel bands seen from 
Copenhagen. 

Faint bands seen low on horizon at 
Newton Stewart. Veil and bands 
seen from Copenhagen. 

No NLC 

Extensive display of faint 
white bands and billows. 


bluish- 


Southern boundary about 55°N. 


Faint NLC seen low on northern 
horizon from Prestwick. 

No NLC 

Cloudy 

Cloudy over British Isles. 
Moderately bright veil and bands 
seen from Denmark through low 
cloud. 

Cloudy 


185 


Station Max. 


: : Limiting 
position elev. 


azimuths 
degrees 


N 1S°E 


$7°N 2°W 
56°S°N 3°W 
56°N 3°W 


55-S°N 1-5°W 


54°N 0-5°W 
$3-S°N 0° 


565°N 3° W 
56°N 30°W 
(aircraft) 
5S6°N 15°W 
(aircraft) 
S6°N 12:S°E 


57-S°N 3-:5°W 


55:5°N 4:5°W 


55-5°N 3°W 
55-S°N45°W 
(aircraft) 

56°N 12:5°E 
55-S°N 3°W 


$7-S°N 3-S°W 


5S°N 4-5°W 


57-5°N 3-5°W 
56:5°N 3°E 
(aircraft) 

56 °N Ww 
55:S°N 4-5°W 


56°N 10°E 





186 Meteorological Magazine, 99, 1970 


NLC were observed on 21 nights during 1969. During the previous two 
summers they were seen on more than 30 nights.*:* Only part of this decrease 
in frequency can be attributed to the greater prevalence of ordinary clouds 
at night during the summer of 1969. The first display occurred a few days 
later and the last, on 20-21 July, at least 10 days earlier than in the previous 
five years. The majority of the displays were faint, only those of 27-28 and 
29-30 June and 20-21 July being of moderate brightness. Whirls were seen 
only during the first two of these three brighter displays. 


We are grateful to the many observers whose reports have been used in 
making this analysis. These synoptic studies are continuing and we invite 
the co-operation of observers who may be prepared to contribute to them. 
Notes on observation of NLC appeared in the Meteorological Magazine, June 
1967, p. 189. An International noctilucent cloud observation manual has just been 
prepared and has been published this year under the auspices of the World 
Meteorological Organization, Geneva. Observations made in western Europe 
should be sent to the Balfour Stewart Laboratory, University of Edinburgh, 
Drummond Street, Edinburgh EH8 gUA, Scotland. 
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REVIEW 


Frederick C. Fuglister. Sixtieth anniversary volume. Deep-Sea Research, Supplement 
to Volume 16. 255 mm x165 mm, pp. v+470, illus., Pergamon Press, 
Headington Hill Hall, Oxford, 1969. Price: £5 (paper-back). 

Thirty-seven papers, mainly on physical oceanography and written with a 

few exceptions by American ‘east coast’ oceanographers, are combined in 

this special volume of Deep-Sea Research to celebrate the sixtieth birthday 
of Frederick Fuglister, the friend and colleague of the contributors and, since 
his conversion from art to oceanography, a dedicated and successful student 
of the Gulf Stream. A surprising proportion of the contributions are of interest 
to meteorologists as well as oceanographers; some are specifically concerned 
with the atmosphere. Haurwitz and Fogle provide a useful review of “Wave 
forms in noctilucent clouds’ (though it lacks reference to some recent work 
on billows); Namias’s paper ‘Autumnal variations in the North Pacific and 
North Atlantic anticyclones as manifestations of air-sea interactions’ 
was written too early to benefit from the recent detailed documentation by 
the Synoptic Climatology Branch; Joanne Simpson’s review ‘On some 
aspects of sea~air interaction in middle latitudes’ emphasizes the role of con- 
vection in connecting the friction layer with synoptic-scale motions. Theoretical 
treatments primarily directed at explaining the structure and movement of 
the ocean frequently have a direct relevance to the atmosphere, for the relevant 
fluid dynamical quantities (the dimensionless numbers named after Rossby, 
Richardson, Brunt—Vaisala and Reynolds) are similar in the two fluids; 
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contributions in this category include Faller and Kaylor on ‘Oscillatory and 
transitory Ekman boundary layers’ (assumed adiabatic and having either 
uniform or a quadratic eddy viscosity), and Geisler and Kraus, who treat 
‘The well-mixed Ekman boundary layer’ as a rigid slab and so neglect the 
turning of the wind with height. 


The papers on experimental (or rather, exploratory, since cruises are usually 
planned with the aim of measuring ocean quantities rather than testing 
hypotheses) oceanography reflect the present transition from the classical 
problem of establishing the general circulation and ocean climate, to the 
growing interest in variability, the oceanographer’s equivalent of weather. 
Of course, there is no simple demarcation between these two aspects of physical 
oceanography and many papers in this volume contain elements of both. 
Especially interesting to meteorologists (since the authors tackle phenomena 
whose atmospheric equivalent has been rather neglected) are the two papers, 
one by Fofonoff the other by Webster, which examine the fluctuations in 
current meter records caused by internal waves (although Webster chooses 
to cover all such variations under the general name ‘turbulence’; better, 
in my opinion, to avoid this confusing definition, which arises from the 
regrettably pervasive —5/3 power law in geophysical spectra). 

Taken together, these varied contributions combine to reveal the richness 
of opportunity offered by physical oceanography, and how profound is our 
ignorance of the ocean. Worthington illustrates the latter with his admirably 
honest paper ‘An attempt to measure the volume transport of Norwegian 
Sea overflow water through the Denmark Strait’. This costly experiment, 
involving a sophisticated array of moored current meters, ended with most 
of the apparatus strewn in an irrecoverable tangle on the sea floor. Happily, 
Worthington was able to reconstruct what had happened from the record 
of one surviving current meter located on the edge of the Strait; apparently 
the overflow is not a steady weak current but occurs as sporadic torrents of 
unprecedented strength, one of which destroyed his apparatus. 


Meteorologists and oceanographers will both benefit from reading the 
many excellent papers in this birthday presentation to Frederick Fuglister. 
J. Woops 


OFFICIAL PUBLICATION 


Geophysical Memoirs 
No. 112 Average temperatures, contour heights and winds at,50 millibars 
over the northern hemisphere. By R. A. Ebdon. 


This memoir uses data mainly for the period 1957-61 to describe the 
climatology of that part of the stratosphere near 50 mb (approximately 
205 km). Charts are included which show average temperatures, contour 
heights and winds and also the variability of each over the northern hemisphere 
during the months of January, April, July and October. Similar charts are 
included for regions north of 45°N for the months of February and March 
in order to assist in the understanding and interpretation of events leading 
up to the ‘final warming’ and the usually rapid breakdown of the winter 





188 Meteorological Magazine, 99, 1970 


régime in high latitudes. The very different thermal régimes and the often 
bimodal character of the temperature distribution in high latitudes in winter 
are described. 


The approximately ‘26-month oscillation’ in tropical stratospheric winds, 
i.e. the alternation between easterly and westerly régimes, is described. The 
elliptical nature of the wind distributions at the 50-mb level over much of 
the hemisphere is discussed. 
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